Inelastic lifetimes of hot electrons in real metals 
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We report a first-principles description of inelastic lifetimes of excited electrons in real Cu and 
Al, which we compute, within the GW approximation of many-body theory, from the knowledge 
of the self-energy of the excited quasiparticle. Our full band-structure calculations indicate that 
actual lifetimes are the result of a delicate balance between localization, density of states, screening, 
and Fermi-surface topology. A major contribution from d-electrons participating in the screening of 
electron-electron interactions yields lifetimes of excited electrons in copper that are larger than those 
of electrons in a free-electron gas with the electron density equal to that of valence {4:S^) electrons. 
In aluminum, a simple metal with no d-bands, splitting of the band structure over the Fermi level 
results in electron lifetimes that are smaller than those of electrons in a free-electron gas. 



Electron dynamics in metals are well-known to play 
an important role in a variety of physical and chemical 
phenomena, and low-energy excited electrons have been 
used as new probes of many-body decay mechanisms and 
chemical reactivity 0]. Recently, the advent of timc- 
resolved two-photon photoeniission (TR-2PPE) Q| has 
made possible to provide direct measurements of the life- 
time of these so-called hot electrons in copper , other 
noble and transition metals | p|, f erromagnetic solids , 
and high Tc superconductors . Also, ballistic electron 
emission spectroscopy (BEES) has shown to be capable 
of determining hot-electron relaxation times in solid ma- 
terials [0. 

An evaluation of the inelastic lifetime of excited elec- 
trons in the vicinity of the Fermi surface was first re- 
ported by Quinn and Ferrell [p^ , within a many-body 
free-electron description of the solid, showing that it is 
inversely proportional to the square of the energy of the 
quasiparticle measured with respect to the Fermi level. 
Since then, several free-electron calculations of electron- 
electron scattering rates have been performed, within the 
random-phase approximation (RPA) |p^ , p^ an d with in- 
clusion of exchange and correlation effects~M,|l^ . Band 
structure effects were discussed by Quinn [|18| and Adler 
[p^ , and statistical approximations were applied by Tung 
et al and by Penn Nevertheless, there was no 
first-principles calculation of hot-electron lifetimes in real 
solids, and further theoretical work was needed for the 
interpretation of existing lifetime measurements and, in 
particular, to quantitatively account for the interplay be- 
tween band-structure and many-body effects on electron 
relaxation processes. 

In this Letter we report results of a full ab initio eval- 
uation of relaxation lifetimes of excited electrons in real 
solids. First, we expand the one-electron Bloch states in a 



plane-wave basis [p2| , and solve the Kohn-Sham equation 
of density-functional theory (DFT) |2|] by invoking the 
local-density approximation (LDA) for exchange and cor- 
relation . The electron-ion interaction is described by 
means of a non-local, norm-conserving ionic pseudopo- 
tential and we use the one-electron Bloch states 

to evaluate the screened Coulomb interaction within a 
well-defined many-body framework, the RPA |^^. We 
finally evaluate the lifetime from the knowledge of the 
imaginary part of the electron self-energy of the excited 
quasiparticle, which we compute within the so-called GW 
approximation of many-body theory ]27[ |. 

Let us consider an inhomogeneous electron system. 
The damping rate of an excited electron in the state 
(/>o(r) with energy E is obtained as (we use atomic units 



throughout, i. e., e — h = 
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where S(r,r';i?) represents the electron self-energy. In 
the so-called GW approximation, only the first term of 
the expansion of the self-energy in the screened interac- 
tion is considered, and after replacing the Green function 
(G) by the zero order approximation (G^), one finds 



Im I](r, r'; £;) = ^ 0} (r')Im W{r, r'; cc>)0/(r), 
/ 



(2) 



where lu = E — Ef represents the energy transfer, the 
sum is extended over a complete set of final states (/)/(r) 
with energy Ef {Ep < Ef < E), Ep is the Fermi energy, 
and W{r,r';uj) is the screened Coulomb interaction: 



W{v,v'-Lo) = j dr"e-i(r,r",c^)w(r"-r'). 



(3) 
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Here, ■i;(r — r') represents the bare Coulomb interaction, 
and e~^(r,r',u;) is the inverse dielectric fimction of the 
solid, which we evaluate within RPA p^ . 

We introduce Fourier expansions appropriate for peri- 
odic crystals, and find 
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where 



i?o/(q) = / d^r<AS(r)e''i--</./(r), 



(5) 



and where eQ^Q,(q, a;) represent Fourier coefficients of 
the inverse dielectric function of the crystal. G and G' 
are reciprocal lattice vectors, and the integration over q 
is extended over the first Brillouin zone (BZ). 

In particular, if couplings of the wave vector q -t- G to 
wave vectors q -I- G' with G 7^ G', i. e., the so-called 
crystalline local-field effects are neglected, one can write: 

^-i^J_Sr f dqV I^O/(q+G)l'lni[£G,G(q,c.)] 
^'yiBZ ^ |q + G|' |eG.G(q,^)|2 ' 

(6) 

If all one-electron Bloch states entering both the ma- 
trix elements -Bo/(q + G) and the dielectric function 
eG,G(q, ti^) are represented by plane waves, then Eq. (|^) 
exactly coincides with the GW formula for the scattering 
rate of excited electrons in a free-electron gas (FEG), as 
obtained by Quinn and Ferrell and by Ritchie JTif . 
In the case of electrons with energy very near the Fermi 
energy {E w Ep) the phase space available for real tran- 
sitions is simply E — Ep, which yields the well-known 
{E — EpY scaling of the scattering rate. In the high- 
density limit {ts « 1) [p^, one finds O 
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The hot-electron decay in real solids depends on both 
the wave vector k and the band index n of the ini- 
tial Bloch state (jioir) = e*""uk,n(r)- We have evalu- 
ated hot-electron lifetimes along various directions of the 
wave- vector ||30), and have found that scattering rates 
of low-energy electrons are strongly directional depen- 
dent. Since measurements of hot-electron lifetimes have 
been reported as a function of energy, we here focus on 
the evaluation of t~^{E), which we obtain by averag- 
ing r^^(k, rt) over all wave vectors and bands lying in 
the irreducible element of the Brillouin zone (IBZ) with 
the same energy. The results presented below have been 
found to be well converged for all hot-electron energies 
under study {E — Ep = 1.0 — 3.5 eV), and they all have 
been performed with inclusion of conduction bands up to 



a maximum energy of ^ 25 eV above the Fermi level. The 
sampling of the BZ required for the evaluation of both 
the dielectric matrix and the hot-electron decay rate of 
Eqs. (^) and (^ has been performed onl6xl6xl6 
Monkhorst-Pack meshes . The sums in Eqs. (|) and 
(^) have been extended over 15 G vectors of the recipro- 
cal lattice, the magnitude of the maximum momentum 
transfer q -I- G being well over the upper limit of 'iqp 
[qp is the Fermi momentum). 

Our ah initio calculation of the average lifetime t[E) 
of hot electrons in Cu, as obtained from Eq. (^) with full 
inclusion of crystalline local-field effects, is presented in 
Fig. 1 by solid circles. The lifetime of hot electrons in 
a FEG with the electron density equal to that of valence 
(4s^) electrons in copper (r^. = 2.67) is exhibited in the 
same figure, by a solid line. Both calculations have been 
carried out within one and the same many-body frame- 
work, the RPA; thus, the ratio between our calculated ah 
initio lifetimes and the corresponding FEG calculations 
unambiguously establishes the impact of the band struc- 
ture of the crystal on the hot-electron decay. Our ah ini- 
tio calculations indicate that the lifetime of hot electrons 
in copper is, within RPA, larger than that of electrons in 
a FEG with 7-3 — 2.67, this enhancement varying from a 
factor of ~ 2.5 near the Fermi level {E — Ep = 1.0 eV) 
to a factor of ~ 1.5 for E — Ep = 3.5 eV. We have also 
performed calculations of the lifetime of hot electrons in 
copper by just keeping the 4s^ Bloch states as valence 
electrons in the pseudopotential generation. The result 
of this calculation nearly coincides with the FEG calcu- 
lations, showing the key role that d-states play in the 
electron-decay mechanism. 

First of all, we focus on the role that both localiza- 
tion and density of states (DOS) available for real ex- 
citations play in the hot-electron lifetime. Hence, we 
neglect crystalline local-field effects and evaluate hot- 
electron lifetimes from Eq. ^ by replacing the elec- 
tron initial and final states in |i?o/(q + G)| by plane 
waves, and the dielectric function in |eG,G(q, t^^)! ^ by 
that of a FEG with Vg = 2.67. The result we obtain, 
with full inclusion of the band structure of the crystal in 
the evaluation of Im [eG,G(q, ^)], is represented in Fig. 
1 by open circles. Since the states just below the Fermi 
level, which are available for real transitions, have a small 
but significant d-component, they are more localized than 
pure sp-states. Hence, their overlap with states over the 
Fermi level is smaller than in the case of free-electron 
states, and localization results in lifetimes of electrons 
with E — Ep < 2 eV that are slightly larger than pre- 
dicted within the FEG model of the metal (solid line). 
At larger energies this band structure calculation pre- 
dicts a lower lifetime than within the FEG model, due to 
opening of the d-band scattering channel dominating the 
DOS with energies from ~ 2.0 eV below the Fermi level. 

While the excitation of d-electrons diminishes the life- 
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time of electrons with energies E—Ep > 2 eV, d-electrons 
also give rise to additional screening, thus increasing the 
lifetime of all electrons above the Fermi level. That this 
is the case is obvious from our band-structure calcula- 
tion exhibited by full triangles in Fig. 1. This calcula- 
tion is the result obtained from Eq. (^) by still replacing 
hot-electron initial and final states in |i3o/(q+G)|^ by 
plane waves (plane- wave calculation), but including the 
full band structure of the crystal in the evaluation of both 
Im [eG,G(q, w)] and |eG,G(q, The effect of virtual 

interband transitions giving rise to additional screening 
is to increase, for the energies under study, the lifetime 
by a factor of « 3, in qualitative agreement with the 
approximate prediction of Quinn . 

Finally, we investigate band structure effects on hot- 
electron energies and wave functions. We have performed 
band-structure calculations of Eq. with and without 
(see also Eq. (^) the inclusion of crystalline local field 
corrections, and we have found that these corrections are 
negligible for E — Ep > 1.5 eV, while for energies very 
near the Fermi level neglection of these corrections re- 
sults in an overestimation of the lifetime of less than 5%. 
Therefore, differences between our full (solid circles) and 
plane-wave (solid triangles) calculations come from the 
sensitivity of hot-electron initial and final states on the 
band structure of the crystal. When the hot-electron en- 
ergy is well above the Fermi level, these states are very 
nearly plane wave states for most of the orientations of 
the wave vector, and the lifetime is well described by 
plane-wave calculations (solid circles and triangles nearly 
coincide for E — Ep > 2.5 eV). However, in the case of 
hot-electron energies near the Fermi level, initial and fi- 
nal states strongly depend on the orientation of the wave 
vector and on the shape of the Fermi surface. While the 
lifetime of hot electrons with the wave vector along the 
necks of the Fermi surface, in the TL direction, is found 
to be longer than the averaged lifetime by up to 80%, 
flattening of the Fermi surface along the TK direction 
is found to increase the hot-electron decay rate by up 
to 10% (see also Ref. [|9|). Since for most orientations 
the Fermi surface is flattened, Fermi surface shape ef- 
fects tend to decrease the average inelastic lifetime. An 
opposite behaviour occurs for hole-states with a strong 
d-character below the Fermi level, localization of these 
states strongly increasing the hole- lifetime | ]30[ |. 

Our band-structure calculation of hot-electron life- 
times in Al, as obtained from Eq. (^, is exhibited in 
the inset of Fig. 1 by solid circles, together with life- 
times of hot electrons in a FEG with = 2.07 (solid 
line). Since aluminum is a simple metal with no d-bands, 

Im — Cg^g' i'i^ ^) ^^1-^ described within a free-electron 
model, and band structure effects only enter through the 
sensitivity of hot-electron initial and final wave functions 
on the band structure of the crystal. Due to splitting of 
the band structure over the Fermi level new decay chan- 



nels are opened, and band structure effects now tend to 
decrease the hot-electron lifetime by a factor that varies 
from ^ 0.65 near the Fermi level {E — Ep — 1 eV) to a 
factor of 0.75 iov E-Ep = 3 eV. 

Scaled lifetimes of hot electrons in Cu, as determined 
from most recent TR-2PPE experiments [^^0], are rep- 
resented in Fig. 2, together with our calculated lifetimes 
of hot electrons in the real crystal (solid circles) and in 
a FEG with = 2.67 either within the fufl RPA (solid 
line) or with use of Eq. (^ (dashed line) . Though there 
are large discrepancies among results obtained in differ- 
ent laboratories, most experiments give lifetimes that 
are considerably longer than predicted within a free- 
electron description of the metal. At E — Ep < 2 eV, 
our calculations are close to lifetimes recently measured 
by Knoesel et al in the very- low energy range At 
larger electron energies, good agreement between our 
band-structure calculations and experiment is obtained 
for Cu(llO), where no band gap exists in the ky = di- 
rection. However, one must be cautious in the interpreta- 
tion of TR-2PPE lifetime measurements, since they may 
be sensitive to electron transport away from the probed 
surface and also to the presence of the hole left behind in 
the photoexcitation process. In the case of injected elec- 
trons in BEES experiments no hole is present. A careful 
analysis of electron-electron mean free paths in these ex- 
periments has shown a {E ~ Ep)~'^ dependence of 
hot-electron lifetimes in the noble metals, with an overall 
enhancement with respect to those predicted by Eq. (|^) 
by a factor of ^ 2, in agreement with our band-structure 
calculations. We note that our calculated lifetimes (solid 
circles) approximately scale as {E — Ep)~'^, as a result 
of two competing effects. As the energy increases hot- 
electron lifetimes in a FEG (solid line) are known to be 
larger than those predicted by Eq. (Q) (dashed line) , 
and this enhancement of the lifetime is nearly compen- 
sated by the reduction, at energies well over the Fermi 
level, of band structure effects. 

In conclusion, we have performed full RPA band- 
structure calculations of hot-electron inelastic lifetimes 
in real solids, and have demonstrated that decay rates of 
excited electrons strongly depend on details of the elec- 
tronic band structure. In the case of Cu, a subtle compe- 
tition between localization, density of states, screening, 
and Fermi-surface topology results in hot-electron life- 
times that are larger than those of electrons in a FEG, 
and good agreement is obtained, for E — Ep > 2 eV, with 
observed lifetimes in Cu(llO). For Al, interband transi- 
tions over the Fermi level yield hot-electron lifetimes that 
are smaller than those of electrons in a FEG. 
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FIG. 1. Hot-electron lifetimes in Cu. Solid circles repre- 
sent our full ab imtio calculation of t{E), as obtained af- 
ter averaging r(k,n) of Eq. (^ over wave vectors and over 
the band structure for each k. Solid and dashed lines repre- 
sent the lifetime of hot electrons in a EEC with = 2.67, 
as obtained within the full RPA (solid line) and from Eq. 
(^) (dashed line). Open circles represent the result obtained 
from Eq. (^ by replacing hot-electron initial and final states 
in |_Bo/(q+ G)|'^ by plane waves and the dielectric function 
in |eG,G(q,ci;)|~^ by that of a EEC with rs = 2.67, but 
with full inclusion of the band structure in the calculation of 
Im [eG,G(q, a;)]. Full triangles represent the result obtained 
from Eq. (^ by replacing hot-electron initial and final states 
in |_Bo/(q+ G)|'^ by plane waves, but with full inclusion of 
the band structure in the evaluation of both Im [eG,G(q, w)] 
and |eG,G(q, <^)|~'^. The inset exhibits hot-electron lifetimes 
in Al. Solid circles represent our ab initio calculations, and 
the solid line represents lifetimes of hot electrons in a EEC 
with rs — 2.07, as obtained within the full RPA. 

FIG. 2. Scaled lifetimes, r x {E — Ef)^ , of hot electrons in 
Cu, as obtained from Eq. (Q) with full inclusion of crystalline 
local- field effects (solid circles). Solid and dashed lines repre- 
sent scaled lifetimes of hot electrons in a EEC with rs = 2.67, 
as obtained within the full RPA (solid line) and from Eq. (M) 
(dashed line). Experimental lifetimes are taken from Ref. M 
(Cu[100]: squares, Cu[110]: inverted triangles, Cu[lll]: open 
circles), from Ref. Q with ujfot ~ 1.63 eV (diamonds), and 
from Ref. M (triangles). 
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Figure 1 
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